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Abstract. A new direct method of obtaining reductions of the Toda equation 
is described. We find a canonical and complete class of all possible reductions 
under certain assumptions. The resulting equations are ordinary differential- 
difference equations, sometimes referred to as delay-differential equations. The 
representative equation of this class is hypothesized to be a new version of one 
of the classical Painleve equations. The Lax pair associated to this equation 
is obtained, also by reduction. 



1. Introduction 
Many papers have been written on the Toda equation 

Ut = u{v — v) 

where u = Un{t) = u{n,t), v = Vn{t) = v{n,t), Ut = du{n,t)/dt, u = u{n + !,<:), 
u = u{n — 1, t), U = v{n + 1, i), w = v{n — 1, t). The famlHar exponential form of 
the Toda equation 

(PQn 



(1.1) 



exp{-((3„ - Qn-i)] - cxp{-((3„+i - Q„)} (1.2) 
ige of variables 

exp{-(Q„ - 0„-i)/2} (1.3) 



can be obtained through the change of variables 

1 

1 dQn-l 



(1.4) 



2 dt 

In either case, the solutions depend on two independent variables (n, t) . Equations 
that involve iterates in one independent variable and derivatives in another are 
referred to as differential-difference equations. Below, we also use this term for 
equations in which iterates and derivatives in the same variable appear. 

The Toda equation is also an example of a "lattice model" . Such lattice models 
appear in many physical settings, ranging from the study of thermalization in metals 
to the study of cellular neural networks and optical lattices. For these applications, 
it is essential to understand the wide variety of solutions possible. Reductions 
provide one method of extending our knowledge of the space of known solutions. 

We focus on reductions to equations involving only one independent variable 

H,^ = K{H,H,H), H-.R^M. (1.5) 
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where iJ,, = dH/di], H ~ H{rj,t), H_ = H{rj,t). Examples of such reductions 
of the Toda equation were obtained through Lie symmetry analysis by Levi and 
Winternitz 

1.1. Background. Symmetry reductions of differential-difference equations have 
been studied by many authors. A comprehensive review can be found in Early 
results on symmetry reductions kept n fixed while allowing u and t to be deformed 
by a method based on the classical approach for differential equations that was 
developed by Sophus Lie. An interesting variation on this standard approach was 
considered in [5] by allowing n to deform continuously in addition to u and t. 
The resulting equations arising as reductions are ordinary differential-difference 
equations of the form (jl.Sp . 

In [5], Quispel et al found an equation of the form (jl.Sp as a reduction of the 
Kac-van Moerbeke or Volterra equation. They showed that their reduced equa- 
tion becomes the classical first Painleve equation [3] in a continuum limit. This 
equation was described as a "delay-differential''^ version of a Painleve equation by 
Grammaticos et al. [2]. Other such equations were proposed by Grammaticos et 
al. [2] as delay Painleve equations by using the criterion of the singularity con- 
finement method. In j4j, other ordinary differential-difference equations were found 
as reductions of the Toda equation. We show that these are contained in our re- 
sults. Moreover, we show that our results arc complete under the assumptions given 
below. 

1.2. Direct Method. In 1989, Clarkson and Kruskal [1] used a "direct approach"to 
find reductions of the Boussinesq equation and found new reductions which were not 
captured by the classical Lie symmetry approach. This direct approach was later 
shown to be related to "non-classical" symmetries. We develop a direct approach 
to finding reductions of differential-difference equations. 

The most general form for a reduction is 

u{n, t) = C/(n, t, i7(ry), G(r,)), v{n, t) = V{n, t, H{ij),G{7^)), = r^n, t), (1.6) 

where H and G form a coupled system of equations of the form p.5|) . For Equations 
(|1.1[) . it turns out to be sufficient to take the ansatz 

u{n,t) ^ a{n,t) + b{n,t) H{r]) (1.7a) 
v{n,t) = c{n,t) + d{n,t)G{ij) (1.7b) 

where rj = r]{n,t). Central to our argument are the following rules (stated for a, b 
and H for conciseness, but they apply also to c, d and G) 

Rule 1: If a{n,t) = ao{n,t) + b{n,t)T{ri), then we can take F = w.l.o.g. by 

substituting H{r]) i—>- H{r]) ~ F(?7). 
Rule 2: if b{n, t) has the form &(n, t) = &o('^, t) ^{v)i then we can take F = 1 w.l.o.g. 

by substituting H{r]) i— > H{r])/T{ri). 
Rule 3: If r/{n, t) is determined by an equation of the form F(?7) = r]o{n, t), where F 

is invertible, then we can take F(77) = ij w.l.o.g. by substituting ij i— > F~^(77). 



We note that this reduced equation contains both retarded H(ri — 1) and advanced terms 
H{ri+1), whilst the usual terminology limits the usage of the term "delay" to equations containing 
only retarded terms such as H(ri — 1). 



DIRECT "DELAY" REDUCTIONS OF THE TODA EQUATION 



3 



Definition 1.1. Given non-zero, differentiable and invertible functions, r(77), we 
refer to the transformations H{rj) H{r]) — r(rj), H{rj) H{rj) /T{rj), and jy ^ 
r~^(r;) as the reduction transformations on (II. 7p . 

Our main resuh is given in Proposition 12 . 1 1 below. 

1.3. Outline of Results. In this paper, we obtain direct reductions of Equations 
p.ip . The details of our direct method are given in Section 2. In Section 3, we 
show that the ansatze (|1.7p in fact represent the general case and can be assumed 
without loss of generality. A continuum limit of the resulting differential-difference 
equations is given in Section 4. In Section 5, we find corresponding reduction of 
the Lax pair for the Toda equation. Finally, a conclusion rounds off the paper. 

2. Direct Reduction of the Toda Equation 

Proposition 2.1. Suppose that the ansatze |_?. Tp hold. Then the only possible non- 
linear second-order reduction of Equation of the form \1.5\) that is unique up 
to reduction transformations of H . G and rj is given by 

-cqH + H,^ = H{G-G) (2.1a) 

Po-coG + G„==2(i?2-H2) (2.1b) 

where the reduction is given by ri{n,t) ~ v{n) + cr(i), ^{n) being an arbitrary func- 
tion, with 

a(t)^[ ilog(cot + ci) + C2 if coT^O 
y ai)t + ai otherwise 
where Cj , j = 0, . . . ,2, oq, ai are constants and the reductions of u, v are given by 
the following two respective cases 
(1) Case Co ^ 0; 

w(n,t) = ± ^- H(r]) (2.3a) 



Co t + ci 
Po 1 

Co (cot + ci) ' ' co^H-ci 



«(n,t) = -^— ^ - + C3+ /, Giif) (2.3b) 



(2) Case c^^^: 

u{n,t) = ±aoH{'q) (2.4a) 

v{n,t) =poalt+pi + aoG{7j) (2.4b) 

with pq , pi being constants. 

Proof. Under the ansatze (|1.7p . the Toda equation becomes 

at + btH + bTjtHr, = {a + b H){c - c + dG - dG) (2.5a) 

Ct + dtG + drjtG^ = 2{a'^ - + 2 ah H - 2abH + - b^H^) (2.5b) 

In the following we indicate generic functions of rj (which are assumed to be differ- 
entiable, non-zero and invertible) by the notation Tj{ri). 

Since we seek non-linear reduced equations of the form (|1.5p . we require the 
terms HG and HG to be present in the reduced equation (|2.5ap . Therefore, we 
require that 

bvt ^liv) = bd (2.6a) 

dT2iri) = d (2.6b) 
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Consider Equation (j2.6a|) . which implies d = r]tTi{T]). However, by Rule 2, this 
implies ri(77) = 1 and, therefore, d = rjt w.l.o.g. 

Using this in Equation (|2.6bp . wc obtain rjl ~ 'r]tT2{i]). However, by taking a 
change of variables 

f]^n{^(n,t)), where = (1^(0) , 

we can take rft = r}t, w.l.o.g. This implies that i]t — for some differentiable 

function (j{t), and therefore, r](n,t) = I'i'n) + a{t), where ^(n) is an arbitrary 
function. 

Now consider the second equation (j2.5bp . Requiring that the terms and 
both remain in the reduced equation, we find that wc must have 

b^^df^trsiv) ^ b^ = {<j'{t))\ w.l.o.g. 

where we have used Rule 2 once again. Therefore, we have b = ± cr'(t). 

Now consider the linear terms in Equation (|2.5ap . In particular, if we require 
that terms linear in G remain in the reduced equation, we must have 

ad = brjtTiirj) ^ a = ±a'{t)T4{-q) = bViit]) ^ a = 0, w.l.o.g. 

by an application of Rule 1. 

If, on the other hand, the linear term in H on the left side of the equation remains 
in the reduced equation, then 

bt^bdr^iv) a" = {a'fT5{v). 

However, since cr only depends on t, while rj also depends on n, this can only 
hold if Fs is identically constant. Letting this constant be — cg, we find cr"(t) = 
-Co (cr'(i))^. We integrate this ODE for a{t) to find 



a{t) 



^ log(co t + ci) + C2 if Co 7^ 
aot + ai if Co = 



where Cj, j = 0, 1, 2 and ao, ai are arbitrary constants. 

Finally, if the linear term in H on the right side of Equation (j2.5ap remains, then 

b(c- c) ^ brjt (Fg - Fg) 

where we have taken the liberty of writing the arbitrary function of t] on the right 
as the difference of another such function. This leads to 

c-c = a'(r6 -Fe) c = ct'FsCt?) + 7 

where ^(t) is an arbitrary function of t. However, since a' = d, we now have 
c = dre{ri) + j{t). By Rule 2, wc can take Fg = and, therefore, c = 7(t) w.l.o.g. 
The reduced equations (|2.5aj|2.5b)) arc now 

-Co H + H,^ = H (G - G) 
^'^^^ coG + G„ = 2(H' -H') 



[a'{t)f 

Since this equation can only contain coefficients that are functions of 77, we must 
have 7' = pq {<j'Y , where po is a constant. That is. 



— .J ca^ (cot+ci) 



+ C3, if Co 7^ 



po ao ^ + Pi if Co = 

where C3 and pi are constants. □ 
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Remark 2.1. Our reduced equations (|2.1|1 form a system of differential difference 
equations, which evolves on a sequence of domains containing points 

5* = I*, Vo,W, ■■ ■}, 

where if ijq ~ i^{n) + a{t)^ then 770 = i^in + 1) + (j{t). In the interior of any domain 
in 71, where the mapping i'{n) v{n + 1) is defined, we get a semi-infinite chain 
of points CP and a corresponding sequence of domains on which these iterates are 
defined. Since ^{n) is an arbitrary function, we have an infinite-dimensional family 
of reductions. 

Remark 2.2. We note that two cases of non-linear reductions were found by Levi 
and Winternitz j4j. In these cases, the new independent variable, called y in their 
paper, is given respectively by (A) y = t exp(— arj) or (B) y = t — an, where a is 
a constant. Case (A), after taking logy as a new variable is the sub-case cq = 1, 
ci = 0, C2 = of the above result. Case (B), is the sub-case cq = 0, oq = 1, ai = 
of the above result. In both cases, i^{n) = —an. 

3. Generalization of Ansatze 

Here we show how the ansatze (|1.7p represent the general case. Consider the 
general reduction 

u{n, t) = c(n, t, H{ri{n, t)), G{r/{n, t))) (3.1a) 

v{n, t) = d{n, t, H{ri{n, t)), G{r){n, t))) (3.1b) 

Under these transformations, the Toda equation becomes 

ct + ch rjt + cg G^iTit ^ c(d - d) (3.2a) 
dt + dH Hr, lit + da G„ ??* = 2 (c^ - c^) (3.2b) 

For the reduced equation to each contain the non-linear terms in H , G, wc require 

d^ dTi{r],H,G) (3.3a) 
c"^ =c^r2{v,H,G) (3.3b) 

Rewriting Fi and F2 in these equations appropriately, wc can sum each to get 

c = e{t) F3(77, H, G), d = /{t) T^ii^, H, G). 

Therefore, 

CHVt = e{t) F4(7;, H, G) vt = e{t) T^ii^, H, G) fit)Ti{r,, H, G) 
CH^g{n,t)r5{v,H,G) 
=^ h{n,t)Teir],H,G) +k{n,t). 

Similarly, we find 

d = r{n, t) Triri, H, G) + s{n, t). 

By defining Fg and F7 to be new variables, if necessary, we reproduce the ansatze 
assumed earlier as (|1.7p . 
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4. Continuum Limits 

To find continuum limits, it is useful to first convert the system of equations 
(|2.ip to a single scalar equation. 

Lemma 4.1. Define rijj) = 2 \og{2 H{rf)). Then r{rj) satisfies 

r,p, = co(r'^ - 2 Co) + (exp(r) - 2 exp(r) + exp{r)) (4.1) 
Proof. The proof is by direct calculation. 

r(ry) = 2 log(2i7) ^ = 2 = 2 (cq + G - G) 

=> r,„ = 2co(G-G) +4(i?' -ij2) -2(i/2 _^2-) 
which gives the desired result upon using = exp(r)/4. □ 

Assume for simplicity that rj ~ vq + V- Then a continuum limit of Equation (j4.1[) 
is obtained by assuming 

r~e'^w{z), z = erj, i^q ~ 1, cq = t'^ K 

as e ^ 0. Then w{z) satisfies 

Wz^+6iy2 + az + /3 = (4.2) 

where a and (3 are integration constants. Note that this equation is a scaled and 
translated version of the classical first Painlevc equation. 

5. Lax pair 

We apply the reduction found in §2 to the Lax pair for the Toda equation given 
by Wadati and Toda [2: 

uij) + utp + V %p — Xif) (5.1a) 
Tpt = uip — utp (5.1b) 



Using the results of Proposition 12.11 for the case co ^ 0, we find from Equation 
(Oal) 

±cr'i7^+±cr'if^ + |(-^cr' + C3) +a'G^i} = X'4> (5.2) 

In order for this to be a reduced equation containing coefhcients that are only 
functions of 77, we require C3 = and need to choose a definite sign in b = ±cr', say 
the positive sign. Furthermore, we define a new spectral variable C, = X/cr' . Then 
we obtain for (j){ri^ C,) = ip{n, t), 

V't = Ct 0c + ct' (/)^ = a' (co C0C + 077)- 
Hence we get the reduced Lax pair 

+ (G- — ) = C-^ (5.3a) 



CO' 

Co C?ic + ^H'^-Kl (5.3b) 

We note that the character of the Lax pair has changed from a spectral problem to 
a monodromy problem, since now derivatives in C also appear in the linear problem. 
By differentiating Equation (|5.3ap in two different ways, once with respect to C and 
once with respect to 77, and using (|5.3b|) to replace , while using (|5.3ap to replace 
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c/), we can show that the compatibiUty conditions for the hnear system (|5.3p are 
precisely Equations (|2.ip . 

Remark 5.1. We note that the classical first Painleve equation has no explicit 
solutions expressible in terms of previously known functions. There is no reason to 
believe that Equations should possess any such explicit solutions. However, 
the reduced equations (|2.ip do have the solution i7 = 0, G = 0, if po = 0. The 
corresponding solution of the Lax pair (|5.3p is = 0. These results under ly the 
solvability of the reduced equation and the reduced Lax pair. Local existence of 
solutions to respective initial value problems can be inferred from the standard 
existence theorems for differential-difference equations. 

6. Concluding Remarks 

We have developed a new direct method of obtaining reductions of differential- 
difference equations and applied this method to the Toda equation (|l.ip .Thc com- 
pleteness of our results in the class of "reduction transformations" (see Definition 
II. ip is shown in two steps. First, in Section 2, we found the most general possible 
reduction subject to the ansatze (|1.7p . Second, in Section 3, we showed that the 
ansatze assumed in fact represent the general case. A continuum limit of the re- 
duced equations is found to be related to the first Painleve equation in Section 4. 
Finally, in Section 5, we showed that the reduced equations inherit a linear problem 
from the Toda equation. 

The reductions we have found through the direct method appear to be closely 
related to those obtained by the Lie-symmetry approach [4]. However, we note that 
in the latter approach, the iterations of the independent variable rj are assumed to 
be the simple ones described hy rj — rj -\- a, or rj ^ a i], where a is a constant. 
No such assumption is made in our approach. The reduced equations (|2.ip allow 
a more general functional iteration of 77 to appear. We also note that neither a 
continuum limit of these reductions or their linear problem has been found before. 

The inheritance of a linear problem through the reduction indicates that Equa- 
tions (|2.ip are integrable. These results suggest that the system p.ip is analogous 
to the well known reductions of completely integrable partial differential equations, 
namely the classical Painleve equations. However, questions remain open on how 
close such an analogy might be. 

The solutions of the Painleve equations have a characteristic complex analytic 
structure in the complex plane of their independent variable. Formal Laurent ex- 
pansions of the solutions of the Toda equation and one of its reductions has been 
discussed in [5]. However, while Equations (|2.ip form a first-degree, second-order 
system of differential equations in 77, they form a second-degree, second-order sys- 
tem when considered as a system of difference equations. This introduces a multi- 
valuedness into the system when it is solved for the highest iterate of H, which 
makes any proof of meromorphicity of the solutions slightly delicate. Moreover, the 
simultaneous presence of iterates of (f> and derivatives of </> in 77 makes the linear 
system (|5.3p a type of monodromy problem that has not been studied before. 
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